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ABSTRACT 
A  compact  heat  exchanger  testing  facility  was  constructed  at  the 
United  States  Naval  Postgraduate  School.   The  testing  facility  is  used 
to  collect  heat  transfer  and  friction  factor  data  by  the  steady  state 
steam  to  air  technique.   Heat  exchangers  with  frontal  dimensions  of  up 
to  12  inches  by  12  inches  can  be  tested.   The  maximum  pressure  drop  which 
can  be  produced  across  the  air  side  of  a  heat  exchanger  is  2.5  psi.   The 
testing  facility  was  designed  to  obtain  data  for  exchanger  air  side 
Reynolds  numbers  from  500  to  10,000.   A  description  of  the  apparatus,  the 
experimental  method  and  equations  and  an  evaluation  are  given. 
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NOMENCLATURE 

English  Letter  Symbols 

2 
Aa      -  Air  side  total  heat  transfer  area,  ft 

2 
Ac      -  Air  side  free  flow  area,  ft 

9 

Ac  -     Air  side  fin  area,  ft 

fa  ' 

2 
Afr     -  Air  side  frontal  area,  ft 

Af  -  Steam  side  fin  area,  ft2 

As  -  Steam  side  total  heat  transfer  area,  ft2 

Awa  -  Air  side  wall  area,  ft2 

A^j  -  Steam  side  wall  area,  ft 

C  -  Coefficient  of  discharge  for  a  metering  orifice 

Cc  -  Jet  contraction  ratio 

C^      -  Coefficient  of  discharge  of  a  fluid  orifice  for  a  given  value 
of   3 

c       -   Specific  heat  of  air  at  constant  pressure,  BTU/(lbm  °F) 

d       -  Diameter  of  fluid  metering  orifice,  in. 

D       -  Air  orifice  duct  diameter,  in. 

D(j      -  Hydraulic  diameter  of  duct  immediately  downstream  of  test 
core,  in. 

ERROR  -  Error  in  heat  balance  between  steam  side  and  air  side,  percent 

f  -  Fanning  friction  factor  in  test  core 

f^  -  Friction  factor  in  duct  immediately  downstream  of  test  core 

fm  -  Friction  factor  of  a  passage  in  test  core 

F  -  Velocity  of  approach  factor  for  the  fluid  metering  orifice 

F.  -  Thermal  expansion  factor  for  the  fluid  metering  orifice 

gc      -  Proportionality  factor  in  Newton's  Second  Law,  g  =  32.174 
lbm  ft/(lbf  sec2) 

G       -  Test  core  air  mass  velocity,  lbsm/(hr  ft2) 
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h  -  Unit  conductance  for  thermal  convection  heat  transfer 
(film  coefficient),  BTU/(hr  ft2  °F) 

hc  -  Enthalpy  of  the  condensate  leaving  the  test  core,  BTU/lbm 

hs  -  Enthalpy  of  steam,  BTU/lbm 

H  -  Humidity  ratio  for  air,  lbm  water/ lbm  dry  air 

k  -  Unit  thermal  conductivity,  BTU/ (hr  ft2  °F/ft) 

Kc  -  Contraction  coefficient 

K(j  -  Velocity  distribution  coefficient 

Ke  -  Expansion  coefficient 

£  -  Fin  length  equal  to  one  half  of  total  fin  height 

l^  -  Length  of  duct  from  downstream  face  of  test  core  to  downstream 
pressure  tap,  in. 

&w  -  Wall  thickness  between  steam  and  air  side  of  test  core,  in. 

L  -  Test  core  flow  length,  ft. 

m  -  V2  h/k«f 

•  -  Mass  rate  of  air  flow,  lbm/hr. 

P  -  Pressure,  in.  H20,  in.  Hg,  lbf/in2,  lbf/ft2 

P'b  -  Barometric  pressure,  in.  Hg 

Q  -  Heat  transfer  rate,  BTU/hr 

R  -  Gas  constant  for  air,  53.3  ft  lbf/(lbm  °R) 

r^  -  Hydraulic  radius,  (A^L/A^) ,  ft,  (4  r^-hydraulic  diameter) 

•  -  Mass  rate  of  steam  flow,  lb  m/hr 
t  -  Temperature,  °F 

T  -  Absolute  temperature,  °R,  °K 

U  -  Unit  overall  thermal  conductance,  BTU/(hr  °F  ft2) 

V  -  Velocity,  ft/sec 

v  -  Specific  volume,  ft-'/lb 
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•       -  Mass  rate  of  condensate  from  test  core,  lbjm/hr 
wc  m 
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x  -  Ratio  of  pressure  upstream  of  the  fluid  metering  orifice  to 
pressure  differential  across  the  orifice 

X  -  Humidity  correction  to  the  specific  heat  of  air 

Xm  -  Humidity  correction  to  the  density  of  air 

Y  -  Net  expansion  factor  for  a  square-edged  metering  orifice 

Greek  Letter  Symbols 

3  -  Ratio  of  fluid  metering  orifice  diameter  to  duct  diameter 

A  -  Denotes  difference 

6  -  Fin  thickness,  in. 

n  -  Fin  temperature  effectiveness 

n0  -  Total  surface  temperature  effectiveness  for  air  side 

t|s  -  Total  surface  temperature  effectiveness  for  steam  side 

(j~  -  Ratio  of  free-flow  to  frontal  area,  A  /AfR 

yU.  -  Dynamic  viscosity,  lbm/ (hr  ft) 

p  -  Density,  lb  ti/ft.3 

P  -  Average  density  of  air  in  test  core  for  isothermal  test, 


ave 


lb  rJft3 
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p.       -  Mean  density  of  air  in  test  core  for  hot  core  test,  lb;m/ft 


Dimensionless  Groupings 

NRe     -  Reynolds  number,  (4r^G^() 

Ngt     -   Stanton  number,  (h/GCp) 

Npr     -  Prandtl  number,  OUCp/k) 

j       -  Colburn-j  (Ngt  Npr2/3) 

NTU     -  Number  of  heat  transfer  units  for  an  exchanger  (A  U/$  c„) 
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Subscripts 

1  -  Air  upstream  of  core 

2  -  Air  downstream  of  core 
o  -  Air  metering  orifice 

a  -  Air  side  of  test  core 

c  -  Core 

f  -  Fin 

f  -  Air  film  in  test  core 

s  -  Steam  side  of  test  core 

so  -  Steam  metering  orifice 

si  -  Steam  in  top  header 

s2  -  Steam  in  bottom  header 

w  -  Wall  between  air  and  steam  sides  of  test  core 

s  -  Square  fin 

t  -  Triangular  fin 

c  -  Circular  tube 

Force  and  Mass  Units 

lb,  -  Denotes  pounds  mass  in  distinction  to 

m  r 

lbf  -  Denoting  pounds  force 
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1.    Introduction 

Steady  state  tests    transferring  heat  from  condensing  steam  to 
flowing  air  are  used  to  determine  experimentally  the  heat  transfer 
characteristics  for  compact  heat  exchangers.   The  friction  factor  is 
determined  from  the  air  pressure  drop  across  the  heat  exchanger.   Steam 
and  air  are  used  as  the  working  fluids  in  the  tests  because  they  are 
readily  available  and  produce  effects  which  are  easily  analyzed.   The 
heat  transfer  characteristics  are  given  by  Colburn-j  heat  transfer 
modulus  evaluated  from  the  steady  state  steam  to  air  test  (referred  to 
here-in  as  hot  core  test).   The  friction  factor  is  evaluated  from  both 
the  hot  core  test  and  the  isothermal  test  (referred  to  as  cold  core  test) 
The  Colburn-j  heat  transfer  modulus  and  friction  factor  are  normally 
presented  graphically  as  functions  of  the  heat  exchanger  Reynolds  number. 
Development  of  these  three  parameters  is  given  in  the  theory  section. 

In  addition  to  determining  the  Colburn-j  and  friction  factor,  a 
heat  balance  is  performed  on  the  exchanger  to  check  the  consistency  of 
the  experimental  data. 

The  testing  facility  is  similar  in  design  to  one  built  at  Stanford 
University  in  1947.   The  Stanford  facility  was  built  under  a  United 
States  Navy  contract  to  evaluate  finned  surfaces  for  use  in  compact  heat 
exchangers  and  is  described  by  Kays  in  reference  [3]. 

Compact  heat  exchangers  which  are  suitable  for  steady  state  steam 
to  air  testing  have  plate-fin  geometry  as  shown  below. 
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The  steady  state  steam  to  air  testing  facility  was  designed  and 
built  to  accommodate  a  heat  exchanger  with  frontal  dimensions  of  up  to 
12  inches  by  12  inches.   The  facility  is  capable  of  producing  air  side 
Reynolds  numbers  from  500  to  10,000.   An  air  side  pressure  drop  of  2.5 
psi  can  be  obtained. 

The  testing  facility  is  described  and  the  experimental  method  and 
equations  for  data  reduction  are  presented.   The  data  reduction  equations 
are  conveniently  arranged  for  programing  on  a  digital  computer.   A  test 
of  a  heat  exchanger  of  known  characteristics  is  presented  as  an  evalua- 
tion of  the  testing  facility.   The  operating  procedures  and  the  computer 
program  for  data  reduction  with  accompanying  instructions  are  given  in 
the  appendix. 


2.   Theory 

Heat  Transfer  Characteristics.   The  heat  transfer  characteristics 
of  a  compact  heat  exchanger  are  best  represented  by  dimensionless  group- 
ings.  This  generalizes  the  results  so  that  they  may  be  used  for  working 
fluids  other  than  steam  and  air.   The  succeeding  development  is  that 
commonly  used  in  the  field;  see  Kays  and  London  [4].   The  Colburn-j  heat 
transfer  modulus  is  selected  for  this  purpose  as  it  combines  heat  transfer 
characteristics  and  fluid  properties.   The  Colburn-j  modulus  is  defined  as: 


J"**?*-    (cV)(^ 


*A 


where  the  unit  conductance  for  convective  heat  transfer,  h,  must  be  deter- 
mined experimentally.   The  use  of  a  steam  to  air  system  is  one  of  the 
most  direct  methods  of  experimentally  determining  h. 

The  unit  conductance  for  convective  heat  transfer,  h,  for  the  air 
side,  air  being  the  fluid  of  interest,  is  determined  from  energy  balance 
considerations  within  the  exchanger.   Low  pressure  slightly  superheated 
steam,  about  6  psig  with  5  degrees  superheat,  enters  the  top  of  the  heat 
exchanger  steam  side.   The  steam  reaches  its  saturation  temperature  in  a 
very  short  distance  after  entering  the  exchanger.   Constant  temperature 
is  maintained  as  the  steam  condenses  for  the  remainder  of  travel.   The 
air  is  heated  continuously  as  it  passes  through  the  exchanger.   The 
sketch  below  illustrates  the  temperature  variations  in  the  working  fluids 
as  a  function  of  the  heat  transfer  area. 
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An  energy  balance  for  the  differential  area  dA  yields  the  equation 
for  the  heat  transfer  rate  as: 

dq        =      ft  cp  dTair     =         U  d  A  AT 

(heat  transfer    (heat  transfer  rate     (heat  transfer  rate 
rate)  to  the  air)  from  steam  to  air) 

Rearrangement  gives: 

dTair  .  U  d  A 


(Tsteam-Tair)       m  c 
For  small  over-all  changes  in  fluid  properties  and  a  relatively  uniform  flow 
distribution  in  the  heat  exchanger,  the  over-all  unit  conductance,  U, 
may  be  assumed  to  be  constant.   Integration  over  the  total  area  then 
yields: 


/  Tsteam  -  Tair  out  _\  =  _  U  A« 
\  Tsteam  -  Tair  in  J      m  c 


£n 

V  Tsteam  -  Tair  in   /      m  r. 

The  dimensionless  grouping  on  the  right  side  of  the  equation  is  defined 

as  the  Number  of  heat  Transfer  Units,  NTU.   Thus  the  equation  may  be 

written  as: 


Hn     [    Tsteam  -  Tair  out  \  _  _  NTU 
\  Tsteam  -  Tair  in  J 


or 


(Tsteam  -  Tair   out  \   = 
Tsteam  -  Tair   in     / 


-NTU 
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The  over-all  unit  conductance  based  upon  the  total  air  side  heat 
transfer  area  must  include  a  factor  to  account  for  the  temperature 
variation  in  the  fin  areas.   The  over-all  temperature  effectiveness,  n» 
may  be  written  as:  a 

n  =  l  -  —  (l  -nf) 

A       r 


where  Af  is  the  fin  area  andTjf  is  the  fin  temperature  effectiveness , 
The  effectiveness  for  a  fin  of  constant  conduction  cross  section  as 
given  by  Kays  and  London  [4]  is: 

tanh  m£ 


where 


Summing  the  resistances  to  heat  flow  then  yields: 


I 
U 

(over-all       (steam  side      (wall  (air  side 

resistance)     resistance)      resistance)     resistance) 


The  convective  coefficient  for  the  air  side,  h,  may  be  solved  for  giving: 
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A  value  of  hg  equal  to  2000  BTU/(hr  ft2  °F)  is  assumed.   A  fifty 
percent  error  in  this  assumption  results  in  approximately  two  percent 
error  in  the  value  of  h.   Accordingly  the  heat  transfer  coefficient  for 
the  air  side,  h,  may  be  evaluated  by  experimentally  measuring  the  fluid 
temperatures  and  the  air  mass  rate  of  flow. 

Friction  Factor.   The  flow  of  a  fluid  through  an  exchanger  results 
in  a  pressure  differential  caused  by  the  flow  contraction  as  the  entrance, 
the  fluid  friction  in  the  exchanger  and  the  flow  expansion  at  the  exchanger 


exit.   These  pressure  effects  are  graphically  represented  below. 
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To  define  a  meaningful  heat  exchanger  friction  factor,  the  entrance 
and  exit  pressure  effects  must  be  elimated  from  the  total  pressure 
differential.   The  friction  factor  will  then  be  based  solely  upon  the 
mean  wall  shear.   The  problem  of  entrance  and  exit  pressure  effects  has 
been  investigated  both  analytically  and  experimentally  by  Kays. [2]. 

The  effect  of  pressure  loss  in  the  downstream  duct  is  included  in 
the  total  pressure  drop  in  order  that  point  2  may  be  located  downstream 
sufficiently  to  assure  complete  pressure  recovery.   The  total  pressure 
differential  from  point  1  to  point  2  is  then: 

Ptotal  =  APa  +  APc  "  APb  +  ^duct 


The  expression  for  the  pressure  effect  due  to  flow  contraction  at 
the  entrance  is  given  by  Kays  [2]  is: 


APa    _    vl  (\~  <rz) 


2gc;. 
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2 

KcVc3 

2  9c 


29c 


and  for  the  flow  expansion  at  the  exit  is: 

APb=_vL(i-  **) 

with  (J~  denoting  the  exchanger  free  flow  to  frontal  area  ratio.   K  ind 
Ke  are  contraction  and  expansion  coefficients  which  are  functions  of 
exchanger  geometry  and  Reynolds  number.   Noting  that  the  velocity  may  be 
related  to  the  mass  velocity  and  that  the  density  change  from  section  one 
to  a  and  b  to  two  will  be  very  small  compared  to  the  total  density  change, 
the  above  equations  yield: 
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AP   -  — 


Wi-cr2  +kc) 


and 


Apb=db-  (' "  r*  ~  *e) 
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The  pressure  drop  in  the  exchanger  may  be  found  by  writing  Bernoulli's 
equation  in  the  differential  form  to  include  the  Fanning  friction  term. 
Using  the  mass  velocity,  G,  the  equation  as  given  by  Kays  [3]  is: 
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dx 
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or 
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By  defining  the  mean  specific  volume  in  the  exchanger  as: 

•L 


vm  = 


ml 


Lo->     e 

and  integrating  from  pi,  0  to  pu,  L,  the  exchanger  pressure  differential 


is   given  by: 


AP  =  ^~ 
c        2gcP<3 
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ft 


The  mean  specific  volume  is  found  from  its  definition  by  replacing  it 


with  its  value  as  given  in  the  equation  of  state  for  a  perfect  gas: 


V 


=  *-  *r 


and  expressing  temperature  and  pressure  as  functions  of  x.   As  the 
pressure  variation  is  small  compared  to  the  temperature  variation,  an 
arithmetic  average  value  of  P  may  be  used: 

P  =  P«  +  p*> 
2 

The  expression  for  temperature  as  a  function  of  x  may  be  written  from 

the  temperature-NTU  relation  as: 


NTU 


T  =  Ts  -  (Ts  -  Tfe)  e 


where  T  is  the  steam  temperature.  Substitution  of  the  expressions  for 
temperature  and  pressure  into  the  equation  for  the  mean  specific  volume 
yields:  y^L 


vm  = 


m 


2R 


0- 


NTU 


Ts-(TS-Ta)e"  L 


dx 


Integration  and  using  the  equation  of  state  for  a  perfect  gas  to  simplify, 
and  assuming  no  temperature  change  from  one  to  a  and  b  to  two  the  follow- 
ing equation  obtains: 


Pm       Tsfe+-TbPb 


T.-   U? 


0>-T|)" 


NTU 


An  alternate  expression  may  be  written  by  substituting  the  NTU-tempera- 
ture  relation  into  the  above  equation.   This  expression  is: 


J j   2R    (ts  —  lmTd) 


with  LMTD  being  the  Log  Mean  Temperature  Difference. 


The  pressure  loss  in  the  downstream  duct  may  be  evaluated  as: 

Apducf   _   4  (fdjkj  \     Vd2 

?d  v  Dd ;  2gc 

where  the  density,   p^,  the  friction  factor,  f^,  the  length  of  duct 
from  section  b  to  2,  SL^,    the  hydraulic  diameter,  D^,  and  the  velocity, 
V^,  are  all  quantities  associated  with  the  downstream  duct.   The  velo- 
city in  the  downstream  duct  may  be  related  to  quantities  appearing  in 
the  previous  pressure  loss  equations.   If  the  downstream  duct  has  the 
same  dimensions  as  the  exchanger  and  small  density  variations  are 
neglected,  the  continuity  equation  gives 
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Thence 
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Again  assuming  that  the  density  changes  between  section  1  to  a  and 


b  to  2  are  negligible,  the  various  pressure  drops  are  summed  yielding: 

z  r(i-<r2+KC)  ^  2  ( e,  -  e 2) 
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P-CrZ-K.)  +4^(T2 


Manipulation  of  the  equation  gives: 


Ap    -  ^ 
^tota.  "  2g"c 


fL 


rhe 


+ 


m 


e2 

Kc-((T2+0 

eT-'- 


+• 


'  l  +  Ke+cr2(l  +4pd  - 


ft 


and  solving  the  equation  for  the  friction  factor  we  obtain: 


f  -    rhfrn 
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ptotai      _    (kc-  crz  -1) 

.G2/2gc  e, 
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Therefore  the  friction  factor,  f,  for  a  heat  exchanger  of  given  geometry 
may  be  found  by  experimentally  determining  the  mass  rate  of  air  flow, 
the  entrance  and  exit  air  pressures  and  temperatures  and  the  steam 
temperature. 

Reynolds  Number.   For  uniformity  of  reporting  test  results,  the 
Reynolds  number  is  defined  in  the  same  manner  as  given  by  Kays  and 
London. [4].   This  definition  embodies  the  hydraulic  diameter  as  the 
characteristic  length  and  the  velocity  based  on  the  minimum  free  flow 
area.   Thus: 

Nn     =     DhVe     =     4rhm      _    4rhG 

To  evaluate  the  foregoing  characteristics,  the  test  apparatus 
must  be  capable  of  accurately  measuring  the  air  state  upstream  and 
downstream  of  the  exchanger,  the  steam  state  within  the  exchanger  and 
the  mass  rate  of  air  flow. 

3.    Description  of  the  Test  Apparatus 

General.   Air  is  induced  through  one  side  of  the  compact  heat 
exchanger  by  a  centrifugal  compressor  and  associated  ducting.   Steam 
is  passed  through  the  other  side  of  the  exchanger  by  the  steam  system. 
The  air  and  steam  systems  are  shown  diagramatically  in  Fig.  4  and  Fig.  5. 
The  steam  system  is  designed  to  give  close  control  of  the  steam  state 


within  the  core.   Pressure  and  temperature  instrumentation  are  pro- 
vided for  measuring  the  air  properties  upstream  and  downstream  of  the 
core,  the  mass  rate  of  air  flow,  the  condition  of  the  steam  entering 
and  leaving  the  core  and  the  mass  rate  of  excess  steam  and  condensate 
leaving  the  core. 

The  test  apparatus  is  shown  pictorially  in  Figs.  1,  2,  and  3. 

Air  System.   The  air  system  is  made  up  of  16  gage  galvanized 
steel  ducts  joined  with  1/4  inch  steel  flanges.   The  ducting  prior  to 
the  test  section  has  a  maximum  dimension  of  12  inches  square  and  the 
ducting  after  the  test  section  is  14  inches  in  diameter.   The  lengths 
of  ducting  were  designed  in  accordance  with  the  ASME  Power  Test  Code 
[7]  for  establishment  of  uniform  flow  and  are  shown  in  Fig.  4. 

The  entrance  to  the  test  duct  is  on  the  exterior  of  the  building 
to  minimize  thermal  gradients  in  the  air  stream.   The  entrance  is  three 
feet  square  and  is  covered  by  a  fine  mesh  screen  to  eliminate  foreign 
matter  from  the  air  stream.   The  transition  from  three  feet  square  to 
12  inches  square  is  made  in  the  entrance  section,  each  side  of  which  has 
the  curvature  of  a  quarter  ellipse.   After  a  five  foot  length  of  duct, 
a  second  transition  is  made  from  12  inches  square  to  the  test  core 
dimensions.   This  is  followed  by  the  forward  instrument  section  which 
contains  four  thermocouple  taps,  a  piezometer  ring  and  two  pitot  tube 
taps  for  conducting  vertical  and  horizontal  velocity  surveys.   Follow- 
ing the  test  core  is  the  after  instrument  section  which  contains  nine 
thermocouple  taps  and  a  piezometer  ring.   A  transition  is  then  made  to 
a  14  inch  outside  diameter  duct  at  the  end  of  which  is  a  standard  ASME 
square  edged  orifice.   Piezometer  rings  are  located  one  diameter 
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upstream  of  the  orifice  and  one-half  diameter  downstream  of  the  orifice. 
The  duct  section  downstream  of  the  orifice  also  contains  a  thermocouple 
tap. 

Air  flow  is  regulated  by  a  manually  operated  double  sliding  plate 
valve  located  at  the  turbo-compressor  inlet.   The  valve  plates  slide  on 
Teflon  runners  and  the  drive  screw  is  ball  bearing  mounted  to  eliminate 
maintenance  problems.   The  double  plates  constrict  the  flow  symmetrically 
and  facilitate  fine  control  of  the  air  flow.   Coarse  control  of  the  air 
flow  is  accomplished  with  the  blast  gate  (discharge  valve)  of  the  cen- 
trifugal compressor. 

Air  is  induced  through  the  test  section  ducting  by  a  centrifugal 
compressor  rated  6000  cfm  at  2.5  psi  back  pressure.   The  compressor  is 
driven  by  a  three  phase,  220  volt,  100  hp  motor.   The  compressor  dis- 
charges air  through  a  20  inch  square  duct  to  the  exterior  of  the 
building. 

The  air  flow  is  metered  by  means  of  a  standard  ASME  square  edged 
orifice  designed  to  the  specifications  given  in  the  ASME  Power  Test 
Code  [7] .   Four  orifice  plates  made  of  1/4  inch  type  304  stainless  steel 
were  designed  with  the  following  orifice  diameters.   The  preceeding  duct 
has  an  inside  diameter,  D,  of  13.875  inches. 

d  (in)  3  =  i. 

10.406  0.75 

6.244  0.45 

3.468  0.25 

2.081  0.15 

These  diameters  were  selected  so  that  flow  rates  of  from  250  lbs  per  hour 

to  22,000  lbs  per  hour  can  be  measured  with  the  ranges  overlapping.   The 

orifice  plates  are  fitted  with  a  centering  grove  on  the  bottom  and  a 
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scribe  mark  on  the  top  so  that  they  may  be  accurately  positioned. 

The  ducting  downstream  of  the  test  section  is  insulated  with  two 
inch  fiber  glass  insulation  covered  by  aluminum  foil,  and  the  air  duct 
is  insulated  from  the  test  core  by  1/8  inch  Teflon  gaskets. 

Steam  System.   Heating  of  the  test  core  is  accomplished  by  passing 
saturated  steam  through  one  side  of  the  core.   Before  entering  the  core 
the  steam  is  reduced  to  a  low  pressure,  about  six  psig,  and  a  slightly 
superheated  state,  about  five  degrees  superheat.   In  a  superheated 
condition,  the  state  of  the  steam  can  accurately  be  established  by 
pressure  and  temperature.   The  few  degrees  of  superheat  are  quickly 
lost  in  the  test  core  and  have  negligible  effect  on  the  calculations. 
A  considerable  excess  of  steam  is  passed  through  the  core  to  prevent 
a  thick  film  of  condensate  from  building  up  on  the  heat  transfer  sur- 
faces. 

Steam  is  supplied  from  the  central  steam  heating  system  at  125 
psig  in  a  saturated  state.   It  is  strained  and  then  reduced  to  about 
30  psig  through  a  1  -  1/4  inch  air  operated,  pilot  controlled  pressure 
reducer.   At  this  point  water  is  injected,  if  necessary,  to  desuperheat 
the  steam  and  then  it  flows  through  a  1  -  1/4  inch  centrifugal  separator 
which  removes  all  entrained  moisture.   The  steam  then  passes  through  a 
second  1  -   1/4  inch  air  operated,  pilot  controlled  pressure  reducer 
which  reduces  it  to  about  six  psig.   This  second  reduction  slightly 
superheats  the  steam.   Just  prior  to  entry  into  the  core  the  steam 
passes  through  a  short  transition  section.   To  provide  uniformity  of 
flow  in  the  test  core,  the  transition  section  is  filled  with  aluminum 
shavings  held  in  place  by  a  fine  mesh  stainless  steel  screen.   Steam  and 
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condensate  from  the  bottom  of  the  test  core  flow  through  a  second 
transition  section  to  a  1  -  1/4  inch  centrifugal  separator  where  the 
condensate  is  removed.   A  float  type  trap  receives  the  condensate  and 
discharges  it  at  a  realitively  constant  flow  rate.   The  cbndensate  is 
then  subcooled  in  a  small  tap  water  counter  flow  heat  exchanger  and 
collected  in  a  bucket  for  weighing.   The  exhaust  steam  flows  from  the 
separator  through  a  standard  ASME  square  edged  orifice  with  flange 
pressure  taps  and  is  discharged  to  the  atmosphere. 

Two  orifice  plates  were  selected  to  measure  the  steam  flow  rates 
and  have  the  following  orifice  diameters.   The  inside  pipe  diameter, 

Ds,  is  equal  to  1.25  inches. 

ds 
ds  (in)        6  =  F- 

us 

0.700  .560 

0.890  .712 

The  metering  orifice  was  installed  as  specified  in  the  ASME  Power 

Test  Code  |7|. 

Pressure  and  thermocouple  taps  were  installed  in  both  the  upper 
and  lower  transition  sections  and  a  thermocouple  tap  was  installed 
upstream  of  the  steam  orifice.   The  entire  steam  system  is  well  insu- 
lated to  minimize  heat  loss. 

Compressed  air  for  actuation  of  the  two  steam  pressure  reducers 
is  supplied  by  two  air  pressure  regulators  located  on  the  instrument 
panel.   This  combination  is  designed  to  hold  the  downstream  steam 
pressure  to  _  0.2  psig.   An  air  line  ftom  the  air  supply  of  the  first 
pressure  reducer  was  fitted  into  the  steam  strainer  clean-out  plug  so 
that  the  entire  steam  system  can  be  blown  dry  upon  completion  of  a 
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series  of  test  runs.   This  should  minimize  rusting  when  the  system  stands 
idle  for  long  periods  of  time. 
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Pressure  Instrumentation.   Pressure  instrumentation  is  provided 
for  measuring   the  gage  pressure  in  the  air  duct  upstream  of  the  test 
core  and  upstream  of  the  air  orifice,  the  air  pressure  differential 
across  the  test  core  and  air  orifice,  the  steam  header  pressures  and 
the  pressure  differential  across  the  steam  orifice.   Air  pressures  are 
measured  with  well-type  single  leg  water  manometers  and  inclined  draft 
gages.   Steam  pressures  are  measured  with  well  type  single  leg  mercury 
manometers. 

A  one  inch  and  three  inch  inclined  draft  gage  and  a  60  inch  water 
manometer  are  used  for  measuring  the  air  pressure  differential  across 
the  test  core.   A  30  inch  water  manometer  is  used  for  measuring  gage 
pressure  upstream  of  the  test  core.   The  pressure  tap  downstream  of  the 
test  core  is  located  so  that  full  pressure  recovery  of  the  air  stream 
is  attained.   The  air  orifice  pressure  differential  is  measured  with  a 
three  inch  inclined  draft  gage  and  a  30  inch  water  manometer.   A  60  inch 
water  manometer  is  used  to  measure  the  gage  pressure  upstream  of  the 
air  orifice.   Each  pressure  tap  on  the  air  duct  consists  of  four  1/16 
inch  holes  drilled  symmetrically  around  the  duct  connected  together  by 
a  piezometer  ring  of  1/4  inch  copper  tubing.   One  quarter  inch  plastic 
tubing  is  used  to  connect  the  piezometer  rings  to  the  manometers  and 
draft  gages  and  each  manometer  and  draft  gage  is  isolated  by  a  brass 
toggle  valve. 

Three  30  inch  mercury  manometers  are  used  to  measure  the  gage 
pressure  of  the  top  and  bottom  steam  headers  and  the  pressure  differential 
across  the  steam  orifice.   Water  is  used  as  the  intermediate  fluid 
between  the  steam  and  mercury  and  a  constant  head  of  water  is  maintained 
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by  using  water  pots.   The  change  of  level  in  the  water  pots  is  the  same 
order  of  magnitude  as  the  change  in  level  in  the  mercury  wells  and  these 
changes  have  insignificant  effect  on  the  pressure  readings.   The  water 
pots  are  stainless  steel  cylinders  three  inches  in  diameter  by  three 
inches  high  and  are  located  above  the  top  of  the  manometers  so  that  air 
pockets  will  not  form  in  the  water  lines  connecting  them  to  the  manometers. 
The  connecting  water  lines  are  of  1/4  inch  stainless  steel  tubing  and 
each  manometer  is  isolated  by  a  stainless  steel  toggle  valve,,   The  water 
pots  are  connected  to  the  steam  headers  and  steam  orifice  by  5/8  inch 
copper  tubing.   The  large  diameter  tubing  was  selected  so  that  conden- 
sate can  flow  freely  back  into  the  steam  system  and  not  form  pockets 
which  would  change  the  static  head  on  the  manometers. 

Temperature  Instrumentation.   All  temperatures  are  measured  with 
copper-constantan  thermocouples.   The  thermocouples  are  in  assemblies 
manufactured  by  Honeywell  under  the  trade  name  MegapaKo   The  MegapaK 
consists  of  a  "head"  for  connepting  the  extension  wires,  a  "sheath"  of 
1/8  inch  stainless  steel  tubing  with  the  insulated  leads  extending  down 
through  the  center  and  the  measuring  junction  at  the  end  of  the  sheath. 
The  thermocouples  used  for  measuring  the  air  stream  temperatures  are 
"exposed,"  meaning  that  the  junction  is  extended  one  sheath  diameter 
past  the  end  of  the  sheath.   Those  used  for  measuring  the  steam  tempera- 
tures are  "remote,"  meaning  that  the  junction  is  one  sheath  diameter 
from  the  end  of  the  sheath  and  the  end  of  the  sheat  is  sealed  for  pressure 
tightness.   All  extension  wire  is  polyvinyl  covered  24  gage  copper-con- 
stantan. 

The  thermocouple  voltages  are  measured  with  a  Honeywell  24  point 
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Multipoint  recorder.   The  reference  junction  is  an  ice  junction  and  is 
included  in  the  measuring  circuit.   The  recorder  is  a  D.C.  millivolt 
recorder  with  a  scale  of  zero  to  ten  millivolts.   The  smallest  division 
on  the  chart  paper  is  1/10  of  a  millivolt  and  temperatures  may  be 
estimated  to  1/100  of  a  millivolt,   The  recorder  has  several  modes  of 
operation.   On  the  Multipoint  setting  it  will  balance  and  print  each  of 
the  24  points  consecutively  at  a  rate  of  about  two  seconds  per  point. 
On  the  Selecto-Print  setting  the  recorder  will  balance  and  print  at 
about  two  second  intervals  any  number  of  the  24  points.   On  the  Hold-on 
setting  it  continuously  balances  on  any  one  of  the  24  points  but  will 
not  print.   This  setting  may  be  used  for  obtaining  a  pen  trace  by 
attaching  a  clip-on  pen  to  the  print  mechanism.   The  recorder   should 
be  adjusted  for  zero  and  span  before  use  and  the  adjustment  procedure 
is  given  in  the  operating  procedures. 

The  temperatures  of  the  air  upstream  and  downstream  of  the  core, 
and  the  downstream  of  the  orifice  are  recorded.   Steam  temperature  in 
the  top  and  bottom  headers  and  upstream  of  the  steam  orifice  are  also 
recorded.   To  obtain  an  average  temperature  of  air  upstream  and  down- 
stream of  the  core  a  system  of  thermocouples  is  used.   Upstream  of  the 
core  this  consists  of  three  thermocouples  placed  on  the  horizontal  center 
line.   Downstream  of  the  core  nine  thermocouples  are  arranged  in  a  three 
by  three  grid  with  equal  spacing  between  each  of  the  thermocouples  and 
the  thermocouples  and  the  wall.   In  this  manner  small  variations  in 
temperature  across  the  stream  can  be  averaged  out.   All  other  tempera- 
tures are  measured  by  a  single  thermocouple.   Each  penetration  through 
the  duct  for  a  thermocouple  is  fitted  with  a  Swagelok  compression  fitting 
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so  that  the  penetration  length  can  be  adjusted. 

The  thermocouples  are  numbered  as  shown  in  Figs.  4  and  5  and  the 
same  number  is  printed  beside  the  temperature  reading  on  the  multipoint 
recorder  chart. 

4.    Experimental  Method  and  Equations 

General.   The  experimental  method  and  equations  used  are  essentially 
those  outlined  by  Kays  [3]  for  a  similar  facility.   These  equations  are 
suited  for  use  with  tables  and  graphs.   The  equations  outlined  herein  are 
suited  for  use  in  a  digital  computer  program  which  will  accept  raw  test 
data  and  produce  finished  results.   The  equations  used  in  preparation  of 
original  graphs  and  tables  were  used  if  available  and  approximate 
equations  were  written  for  those  not  availably.   Graphs  of  the  approxi- 
mations are  shown  in  the  appendix  and  test  data  should  be  checked  for 
applicability  of  the  approximations. 

Mass  Rate  of  Air  Flow.   The  mass  rate  of  air  flow,  m,  through  an 
ASME  square  edged  orifice  is  given  in  the  Power  Test  Code  [7]  as: 

m  =359  C  Fd2FAYV4P0p0  '        lbm/hr 

where   C     =  coefficient  of  discharge 

F     =  velocity  of  approach  factor 

d     =  orifice  diameter,  in. 

Fa    =  orifice  thermal  expansion  factor 

Y     =  net  expansion  facto): 

=  pressure  differential  across  the  orifice,  in.  H20 


&o 


■  specific  weight  of  the  air  upstream  of  the  orifice, 
lbm/ft3 
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The  ratio  of  the  orifice  diameter  to  the  duct  diameter  is 

3   =  i 

D 
where        D  =  duct  diameter 

=  13.875  in. 

and  the  velocity  of  approach  factor  is: 


F   = 


yr  <3' 


The  absolute  duct  pressure  upstream  of  the  orifice  is  given  by: 

Po  =  Pb  CO,  490)  -  P0  (0.0361)   psid 

where        P^  =  barometric  pressure,  in.  Hg 

P   =  duct  pressure  upstream  of  the  orifice,  in.  1^0 
(duct  pressure  is  below  atmospheric  pressure; 

A  dimensionless  ratio  of  orifice  pressure  drop  to  orifice  pressure 

is:  APo(0.036l) 

x  =   a-s 

and   the  net   expansion  factor   is    then: 


Y     =      1.0 


_  /JLVq.41  4-  0.3  5  (34) 


The  approximate  equation  for  the  orifice  thermal  expansion  factor 
for  type  304  stainless  steel  gives: 


fa  =  l.o  +(t0-  68)0'S5*IO  5) 


The  coefficient  of  discharge  is  an  unknown  in  the  equation  but  may 
be  quickly  iterated  using  an  equation  suggested  by  Murdock  [6].   The 
equation  is  dependant  upon  the  Reynolds  number  in  the  orifice,  NR  j,  and  is: 

'  I04  \  J/2 


c    -    c0  -+  AC 


Nfted) 
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The  coefficients  C0  and  AC  for  applicable   $'s  from  Murdock  [6]  and  a 
suggested  first  iteration  C  are: 


e 

Co 

AC 

C 

.15 

0.59446 

0.00945 

0.5975 

.25 

0.59483 

0.01037 

0.5966 

.45 

0.59863 

0.01543 

0.6014 

.75 

0.60480 

0.05448 

0.6128 

The  density  of  the  air  upstream  of  the  orifice  may  b,e  found  from 
the  equation  of  state  of  a  perfect  gas  modified  by  a  humidity  correction 
factor  which  gives: 

144  PQXm 
Ko        53.3  4  T0 

where    PQ  =  duct  pressure  upstream  of  the  orifice,  psia 

TQ  =   temperature  of  air  upstream  of  the  orifice,  °R 

^  =  humidity  correction  factor  for  density  as  given  by 
Kays  and  London  [4] 


0  +H^ 


(I  +  1.607  H") 

and      H   =  humidity  ratio,  lbsm  R^O/lb  dry  air 

The  dynamic  viscosity ,  /A.Q,    for  air  at  one  atmosphere  pressure  is 

given  in  reference  (1)  as: 

3/2 

..    _      0.00  3  527   Tok 

Mo  -         -r 

7  Tok  +    110.4 

where    T  ,  =   temperature  of  air  in  duct,  °K 
The  orifice  Reynolds  number  may  be  found  from: 

_  15,2  8  m 
NRed"  /iQ  d 
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Heat  Transfer  Calculations.   The  unit  conductance  for  thermal  con- 
vection heat  transfer  as  given  by  equation  one  in  the  theory  section  is: 

1  


h*    no 


jtittztz  -T-S- 


U  aP1sAshs  !2Awak 


wa   w. 


where     ha   =  unit  conductance  for  thermal  convection  heat 

transfer  to  the  air,  BTU/(hr  ft2  °F) 

y\Q  =  total  surface  effectiveness  for  the  air  side 

U  =  overall  unit  conductance,  BTU/ (hr  ft2  °F) 

Aa  =  total  heat  transfer  area  on  the  air  side,  ft2 

Tls  =  total  surface  effectiveness  for  the  steam  side 

As  =  total  heat  transfer  area  on  the  steam  side,  ft 

h    =  unit  conduction  for  thermal  convection  on  the 
steam  side,  BTU/ (hr  ft2  °F) 

0    =   thickness  of  the  wall  separting  the  steam  side 
from  the  air  side,  in. 

2 

Awa  =   area  of  the  wall  on  the  air  side,  ft 

k^       =      thermal  conductivity  of  the  wall,  BTU/(hr  ft2  °F/ft) 
The  film  coefficient  for  convection  on  the  steam  side,  h  ,  is  unknown 
and  must  be  estimated.   Any  reasonable  estimate  is  considered  adequate 
since  a  50  percent  error  in  the  estimate  results  in  about  two  percent 
error  in  ha.   A  value  of  hs  equal  to  2000  BTU/(hr  ft2  °F)  has  been  used. 

The  specific  heat  of  dry  air  may  be  considered  constant  and  is  assumed 

t0be:  cpdd=  0.241        eTU/(lbm,,F) 

The  humidity  correction  factor  for  the  specific  heat  as  given  by  Kays 
and  London  [4]  is: 


X  = 


+  I.9I5H 


'c  I  +  H 

where  H     =      humidity   ratio,    lb        H.20/lbm  dry   air 
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and  the  corrected  specific  heat  is  then: 

cp  =  cpda  Xc       BTy/(lbm°F) 

The  bulk  average  temperature  in  the  test  core  is: 

tb=    '  I      ^  F 

where     t^  -     air  temperature  upstream  of  the  test  core,  °F 

t2  =  air  temperature  downstream  of  the  test  core,  °F 
and  the  film  average  temperature  is  then: 

t.-    *j£±±l      °f 

r      2 

where     ts  =   the  saturation  temperature  corresponding  to  the 

arithmetic  average  steam  pressure  in  the  core,  °F 

The  number  of  heat  Transfer  Units  may  be  found  from  the  temperature 

NTU  relation  as  given  in  the  theory  section,  yielding: 

-  ■  -£4) 

and  the  overall  .unit  conductance  is  then: 

_    cp4,NTU        BTU/(hr  ft2  -F) 
U  Ad 

For  most  plate-fin  type  geometries  the  equation  for  the  temperature 
effectiveness  of  a  straight  fin  with  constant  conduction  cross  section 
may  be  used  to  a  good  approximation. [3] .   The  fin  temperature  effective- 
ness for  the  steam  side,  rif s »  is: 

tanh(Vnsls/l2y 
%f*   ~        (ms5s/l2^ 

where     &s     =     effective  fin  length  for  the  steam  side  and  is 
equal  to  one-half  the  total  fin  height,  in. 


/  24  hs  '   „  • 
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and       hg   =   2000  BTU/(hr  ft2  °F)  as  previously  estimated 

kfs   =   thermal  conductance  for  the  fin  on  the  steam 
side,  BTU/(hr  ft2  °F/ft). 

6S   =  fin  thickness  on  the  steam  side,  in. 

The  overall  temperature  effectiveness  for  the  steam  side,  ns,  is  then: 

*»•  --  As 

2 
where     A^   =  wall  area  on  the  steam  side,  ft 

2 
Af   =  fin  area  on  the  steam  side,  ft 

2 
As   =   total  heat  transfer  area  on  the  steam  side,  ft 

The  fin  temperature  effectiveness  for  the  air  side,7|ra,  is: 

tanh(m3)la/l2) 

where     £a   =   effective  fin  length  for  the  air  side  and  is  equal 

to  one-half  the  total  fin  height,  in. 


m^  = 


24  h*  '   ^-1 


*-   ft 


a  "  1/  Kfa  ^a 

and       kjs   =   thermal  conductance  for  the  fin  on  the  air  side, 

BTU/(hr  ft2  °F/ft). 

6a   =  fin  thickness  on  the  air  side,  in. 

ha   =  film  coefficient  for  convection  on  the  air  side. 

It  may  be  taken  equal  to  U  for  the  first  iteration 
BTU/(hr  ft2  °F). 

The  overall  temperature  effectiveness  for  the  air  side,  n,  is  given  by: 

AWa  +TlfaAfdi 
^°  Aa 


where     A  =  wall  area  on  the  air  side,  ft 

wa  ' 

A/>a  =  fin  area  on  the  air  side,  ft 

Aa  =  total  heat  transfer  area  on  the  air  side,  ft2 
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The  air  mass  velocity  in  the  core,  G,  is  given  by: 


G  =  A   »bm/(hr  ft2) 


where     A  =   the  free  flow  area  of  the  core  on  the  air  side,  ft'1 
Stanton  number  in  the  core  can  now  be  evaluated  as: 

The  dynamic  viscosity  evaluated  at  the  film  average  temperature  and 
one  atmosphere  pressure  as  given  by  Hilsenrath  [1]  is: 

iA     _    0.003  527    'f.k         Ib^/fhr   ft) 
Mf  '     T*„  +  110-4  m/l  ' 


'f  k 


where 


Tf,  =  film  average  temperature,  °K 


The  dynamic  viscosity  evaluated  at  the  bulk  average  temperature  is: 

3/2 
0.003527    Tt/        lb/(hrft) 

^b         Tbk     +  ll0*4 

where     T,  ,  =  bulk  average  temperature,  °K 

Ther  thermal  conductivity  of  air  at  the  film  average  temperature 
and  one  atmosphere  pressure  as  set  forth  by  Hilsenrath  [1]  is: 

,l/2 
_   0.001529     Tfk BTU 

24  5.4  1  hr  ft*  Vft 


kf  = 


+ 


TfkXI0 


k7t7i 


Where     T^  =   film  average  temperature,  °K 

The  Reynolds  numbers  evaluated  at  film  average  temperature  and  bulk 
average  temperature  are  evaluated  by: 


N 


and 


4fhG 
Ref  JJf 


N 


Reb    ~    yU 


4rbG 
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where  hydraulic  radius  r^  is  defined  by  Kays  and  London  [3]  to  be  based 
on  the  free  flow  area  and  takes  the  form 

ACL 


rh  = 


ft 


where     L  =   core  flow  length,  ft 
Prantdl  number  in  the  core  is: 


N 


Pr 


_  Mr  cp 


where  the  fluid  properties  are  evaluated  at  the   film  average  temperature 
in  the  test  core. 

The  Colburn-j  heat  transfer  modulus  is  then: 

J  E  N^Np/ 
which  concludes  the  heat  transfer  calculations. 

Friction  Factor  Calculations.   The  equation  for  the  core  friction 
factor  as  developed  in  equation  two  in  the  theory  section  may  be  written 


as 


f     = 


rhpm      \    APC (4.3255  X  IP9) 


L 


Kc  -0  +  <J~2) 


K(,  +  i  -t-  <r 


^4     Dd 


where 


Pm 
L 

APC 
G 


=   core  hydraulic  radius,  ft 

mean  density  in  the  core  lb  /ft^ 

J  m 

core  flow  length,  ft 

pressure  drop  across  core,  in.  H?0 

core  mass  velocity,  lbm/hr  ft^ 
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Kc   =  contraction  coefficient 


0~~  =  the  ratio  of  free  flow  area  to  frontal_area 

\\  =  density  of  the  air  upstream  of  core,  lb/ft^ 

Kg  =  expansion  coefficient 

\2  =  density  of  the  air  downstream  of  core,  lbm/ft 

fd   =   friction  factor  for  the  duct  downstream  of  the 
core 

Ift       =      length  of  duct  from  the  core  to  the  downstream 
pressure  tap,  in. 

Del   =  hydraulic  diameter  of  the  duct  downstream  of 
core,  in. 

The  duct  friction  factor,  fd,  is  assumed  constant  at  .015,  a  value 
corresponding,  to  a  nearly  smooth  pipe  for  a  Reynolds  number  range  of 
100,000  to  500,000  as  may  be  expected  in  the  downstream  duct. 

The  mean  density  in  the  case  of  heating  is  given  by: 


em 


T,e,  +T?e 


1M 


2V2 


I 


t  Tl  Z  T2 


NTU 


bmA- 


where 


l2 
Ts 


upstream  temperature,  °R 

=  downstream  temperature,  °R 

=  Saturation  temperature  corresponding  to  the  average 
steam  pressure  in  the  core  °R 


p^   =   upstream  density 

m      0.097318    XmP,  ,bm/ft3 

Tl 

Pi    =  upstream  pressure,  in.  H2O  absolute 

Xm   =  humidity  correction  factor  as  previously  given 

P2   =  downstream  density 

=    0.09318    Xm   Pg  lbm/ft3 

T2 
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•2   =  downstream  pressure,  in.  1^0  absolute 
In  the  isothermal  case,  the  mean  density  is  equal  to  the  arithmetic 
average  density  and  since  the  temperature  remains  constant 

0.097318    Xm(P,trj)     |b„/ft3 
Km  ~   C^ve  2T1 


m 


The  equations  for  the  contraction  coefficient,  Kc,  and  the  expansion 
coefficient,  Ke,  are  developed  by  Kays  [2]  and  are  given  by 


Kr  = 


I  -  2CC+  4  fc  Kjj) 


and 


Ke  =     I  -  2Kd  CT  -+-  (Tl 


where     K^   =  velocity-distribution  coefficient 


<r  -  Ac/a 


fr 


Cc   =  jet  contraction  ratio 
The  jet  contraction  ratio  is  approximated  as 

Cc  =  0.611  -hO.45  cr  +  0.344  <r5'7 

The  velocity-distribution  coefficient  K^  is  dependent  upon  both  core 
geometry  and  film  Reynolds,  NR  f.   The  usual  geometries  used  in  plate 
fin  type  construction  are: 


square 


triangular 


For  laminar  flow,  l%ef  less  than  2000,  Kays  [2]  gives  the  values 


square 


K. 


=  K, 


=   1.39 


triangular    K^   =  K^   -   1.43 
and  for  turbulent  flow 
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square    Kd   =  Kd   =1+1.17  (Kd  -  1) 
triangular    Kd   =  Kd   =   1+1.29  (Kd  -  1) 

where  the  velocity-distribution  coefficient  for  circular  tubes,  Kd  ,  is 

given  by:  

Kdc  =   1.09068  (4fm)   -h  0.0584 -^4fm+  I  ' 

and  _  q    p 

fm  =     0.049  NRef    u'd 

If  the  boundary  layer  is  continually  tripped,  as  in  the  case  of  louvered 
fin  or  off-set  fin  construction,  film  Reynolds  number  should  be  taken 
equal  to  infinity  which  yields  Kd  =  1.0  in  all  cases. 

Heat  Balance  Calculations.   The  heat  transfer  rate  to  the  air, 
Q  air,  is: 

Q  air  =  m  cp  (t2  -  ti)  BTU/hr 
The  total  heat  given  up  by  the  steam,  Q  steam,  is  equal  to  the  heat 
given  up  by  the  excess  steam  as  it  passes  through  the  core  plus  the  heat 
given  up  by  the  condensate  in  the  core. 

Q  steam  =  s  (hi  -  h2g)  +  wc  (hx  -  h2f)  BTU/hr 
where     s  =  mass  rate  of  flow  of  excess  steam,  lbm/hr 

wc  =  mass  rate  of  flow  of  the  condensate,  lbm/hr 

h-L  =  enthalpy  of  the  steam  at  the  top  header  pressure  and 

temperature  BTU/lb 

r  m 

h2  =   enthalpy  of  saturated  steam  at  the  bottom  header 
pressure,  BTU/lbm 

h2f=  enthalpy  of  saturated  liquid  at  the  bottom  header 


pressure,  BTU/lb 


m 


The  mass  rate  of  flow  of  the  excess  steam,  s,  as  measured  by  a  square 
edged,  standard  ASME  orifice  with  flange  pressure  taps  is  found  from  the 


33 


applicable  equation  given  in  the  ASME  Power  Test  Code  [7}: 


S  =  359  Cs  Fs  df  FAsYs  J]^       \  bm/  hr 

where     Cg   =   coefficient  of  discharge 

Fs   =  velocity  of  approach  factor 

d    =  orifice  diameter,  in. . 
s 

^As  =  orifice  thermal  expansion  factor 

Yg   =  net  expansion  factor 

Z\PS   =  pressure  differential  across  the  steam  orifice, 
in.  H20 

Vs   =   specific  volume  of  the  steam  passing  through  the 
orifice,  ft3/lbm 

The  velocity  of  approach  factor,  Fs,  is: 

r-      ,   I 


where 


and 

=   1.25  in. 
The  orifice  thermal  expansion  factor  is  taken  as  constant  corres- 
ponding to  a  steam  temperature  of  220  °F.   For  type  304  stainless  steel 


FAs  =   1.003 
The  net  expansion  factor  for  steam,  Ys,  is: 

YS=|-(0.4I  +0.l5e>?)(*s/\'3) 

where  x  x 

xs=  APs/Ps 

The  upstream  steam  pressure,  P  ,  is  taken  as  the  saturation  pressure 
corresponding  to  the  upstream  temperature  and  is  given  by  the  approximate 


's 

VT^ 

ft4 

Ps 

=    ds 

Ds 

=      steam  pipe 

diameter 

34 


equation 

Ps=  16.5  -f  (tso-2l8X0.376)    psta 

The  specific  volume  is  given  by  the  approximate  equation: 

vs  =  24,OI-O.4(tso-2I8)  ft3/\b 

The  coefficient  of  discharge,  C  ,  is  iterated  by  the  same  method 
that  was  used  for  the  air  orifice.   The  values  of  C  and  AC  for  applicable 
values  of  3q  are: 


6S  ds  C0  =  Cos       AC  =  ACs 

.56  .700  in.  .6085  .025 

.71  .890  in.  .6090  .041 

and  then 

cs=cos+/Wj^_f2 

VNResd> 

The  steam  orifice  Reynolds  number,  Nnescj,  is: 

_    0.004244    S 
NResd         /U$  ds 

and  the  dynamic  viscosity  is  given  by  the  approximate  equation: 

ib 


^5  =  8.25  X  l0"6+(ts0-200)(l.4  5  X  IC-8)  — 


m 


sec 


The  enthalpy  of  the  superheated  steam  in  the  top  header  is  approxi- 
mated by  two  equations,  the  first  giving  the  variation  with  temperature 
at  a  constant  pressure,  Pg4  =  18.0  psia  and  the  second  equation  the 
variation  with  pressure.   These  equations  are: 

hgl  =   1158.0  +  (ts|  -  230.0)  (.5)  BTU/lbm 
where     h  ,   =   enthalpy  at  18.0  psia* 

tsl  =   toP  header  steam  temperatrue,  °F 
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thence    hgl   =  h^  -  (Psl  -  18.0)  (.35)  BTU/lbm 
where     Psi   =   steam  pressure  in  top  header,  psia 
The  enthalpy  of  the  saturated  steam  in  the  bottom  header  is  found 
from  the  approximate  equation: 

hs2   =   1155.3  +  (Ps2  -  19.0)  (.934)  BTU/lbm 
where     P  o   =   stream  pressure  in  bottom  header,  psia 
and  the  enthalpy  of  the  condensate  in  the  bottom  header  is: 
hc    =   193.4  +  (Ps2  -  19.0)  (2.64)  BTU/lbm 
The  mass  rate  of  flow  of  the  condensate,  wc,  is  found  by  timing 
the  condensate  collected. 

5.    Evaluation 

A  steady  state  steam  to  air  test  was  conducted  on  a  compact  heat 
exchanger  manufactured  by  Hamilton  Standard  Corporation.   The  test  core 
had  the  approximate  dimensions  of  a  six  inch  cube.   The  construction  was 
the  plane  plate-fin  type  with  flow  passages  between  fins  being  approxi- 
mately square.   One  side  of  the  core  was  "straight,"  meaning  that  the  flow 
passage  was  straight  for  the  length  of  the  exchanger.   The  other  side  was 
"ruffled,"  meaning  that  the  flow  passage  was  approximately  a  cosine  curve. 
The  plate  spacing  was  0.156  inches  and  the  fin  spacing  was  18  fins  per 
inch  for  both  sides.   All  plate  and  fin  material  was  0.004  inch  stainless 
steel.   A  set  of  "standard"  curves  for  Colburn-j  and  friction  factor  vs. 
Reynolds  number  were  supplied  with  the  test  core.   The  "standard"  curves 
were  for  plate  and  fin  spacing  similar  to  the  core  tested. 

The  test  was  performed  using  the  straight  side  of  the  core  for  air 
and  the  ruffled  side  for  steam.   The  various  components  of  the  facility, 
the  air  system,  steam  system  and  instrumentation  were  checked  out  and 
proved  to  be  very  satisfactory.   There  was  a  problem  in  defining  flow  and 
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heat  transfer  areas  for  the  exchanger  because  structural  stiffeners  had 
been  placed  at  the  end  of  each  plate-fin  passage.   The  stiffeners  were 
placed  at  an  angle  across  the  flow  path  and  thereby  restricted  the 
flow  in  the  peripheral  flow  passages.  This  restriction  would  tend  to 
increase  with  increasing  flow  rates.   The  test  results  based  upon  the 
minimum  free  flow  area,  which  neglects  the  peripheral  flow  passages 
affected  by  the  stiffeners,  was  plotted  on  the  standard  curves  as 
shown  in  Fig.  6. 

The  next  test  should  be  on  a  heat  exchanger  core  with  previously 
determined  characteristics  and  with  given  or  definable  geometry. 
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Figure  6,     Colburn-j  and  friction  factor 
vs«  test  core  Reynolds  number. 
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APPENDIX  I 
OPERATING  PROCEDURE 

General.   The  preliminary  check  of  instrumentation  for  both  the 
hot  core  and  cold  core  tests  is  the  same. 

The  multipoint  recorder  should  be  allowed  to  warm  up  two  hours 
prior  to  calibration.   To  calibrate  the  instrument,  set  the  print 
mechanism  to  "Hold  On"  indicating  any  one  of  the  points  not  used  for  a 
thermocouple.   This  open  circuits  the  thermocouple  measuring  circuit. 
Place  a  D.C.  millivolt  source  across  the  terminals  at  the  back  of  the 
recorder  and  set  the  source  to  one  millivolt.   Adjust  the  recorder  zero 
so  that  one  millivolt  is  indicated.   Set  the  source  to  nine  millivolts 
and  adjust  the  span  so  that  nine  millivolts  are  indicated.   Repeat  the 
zero  and  span  adjustments  until  the  recorder  is  calibrated. 

All  manometers  should  be  adjusted  for  zero  with  the  isolation  valves 
open  prior  to  lighting  off  the  air  or  steam  systems.   The  two  steam 
header  mercury  manometers  must  be  zeroed  at  about  three  inches  due  to  the 
head  of  water  between  the  manometers  and  the  water  pots.   The  water  lines 
between  the  mercury  manometers  and  the  water  pots  are  susceptible  to  the 
formation  of  air  pockets  which  will  cause  erroneous  readings. 
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Cold  Core  Test.   In  order  to  assure  an  even  distribution  of  data 
points,  the  approximate  pressure  differential  across  the  air  orifice 
should  be  calculated  for  each  run.   Calculate  the  mass  rate  of  air 
flow  at  a  core  Reynolds  number  of  500  and  10,000  and  plot  these  two 
points  on  a  log  log  plot.   Divide  the  straight  line  connecting  the  points 
into  equal  increments,  the  number  depending  upon  the  desired  number  of 
runs.   The  orifice  pressure  differential  may  then  be  read  from  Fig.  7 
for  each  mass  rate  of  air  flow. 

Open  all  manometer  valves  and  adjust  for  zero.   Close  all  manometer 
valves  prior  to  lighting  off.   Set  the  multipoint  recorder  to  print  the 
thermocouples  in  use  by  raising  the  corresponding  buttons  on  "Selecto- 
Print."   Turn  the  chart  motor  on  and  select  the  desired  speed.   Close 
the  compressor  blast  gate  (discharge  valve)  and  energize  the  compressor 
motor.   Open  the  blast  gate  to  desired  setting.   Adjust  the  sliding 
plate  valve  for  the  desired  air  orifice  pressure  differential  and  observe 
instrumentation  until  steady  state  obtained. 

Record  the  following  data  for  each  run: 

1.  Atmospheric  pressure  (0.001  in,  Hg) 

2.  Atmospheric  dry  bulb  temperature  (0.1  °F) 

3.  Atmospheric  wet  bulb  temperature  (0.1  F) 

4.  Air  orifice  size 

5.  Orifice  upstream  pressure  (0.01  in,  HL0) 

6.  Orifice  pressure  differential  (0.001  in.  H20  for  APo  <3",  and 
0.01  in.  H20  for  AP&  >3") 

7.  Core  upstream  pressure  (0.01  in.  H20) 

8.  Core  pressure  differential  (0.001  in,  H20  for  APc  <3",  and  0.01 
in.  H20  for  APc  >3"). 
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9.   Orifice  air  temperature  (0.01  millivolt) 
10.   Core  upstream  air  temperature  (0.01  millivolt) 

Items  five  through  ten  are  recorded  three  times  in  the  same  sequence 
and  averaged  on  the  data  sheet. 

Hot  Core  Test.   The  instrumentation  checks  and  the  procedures  for 
lighting  off  the  air  system  are  the  same  as  for  the  Cold  Core  Test. 

The  steam  system  is  activated  by  the  following  procedures.   Obtain 
100  psig  low  pressure  air  to  the  reducer  central  valves.   Turn  on  the 
tap  water  supply  for  the  condensate  cooler  and  the  desuperheater .   Crack 
the  desuperheater  supply  valve.   Open  the  by-pass  valves  for  the  steam 
reducers  and  crack  the  steam  stop  valve.   Allow  the  system  to  warm  up 
slowly.   When  warm-up  is  completed  close  the  by-pass  valves  and  open  the 
steam  stop.   Charge  the  first  steam  reducer  to  about  30  psig  and  the 
second  to  about  six  psig.   Obtain  a  steam  state  in  the  top  steam  header 
of  about  six  psig  with  five  to  ten  degrees  superheat  by  adjusting  the 
first  steam  reducer  and  the  desuperheater.   The  excess  steam  flow  rate 
should  be  at  least  five  times  the  condensate  flow  rate.   Two  steam 
orifice  plates  are  provided;  the  larger  to  be  used  for  high  air  flow  rates. 
Condensate  should  be  collected  for  at  least  five  minutes  at  high  air  flow 
rates  and  for  at  least  ten  minutes  at  low  air  flow  rates. 

Steady  state  must  be  achieved  before  each  run.   This  will  require 
about  one  hour  prior  to  the  first  run  and  15  minutes  between  runs.   Steady 
downstream  core  temperatures  are  the  best  indication  of  steady  state. 

The  following  data  is  recorded  for  each  run: 

1.  Atmospheric  pressure  (0.001  in.  Hg) 

2.  Atmospheric  dry  bulb  temperature  (0.1  °F) 


42 


3.  Atmospheric  wet  bulb  temperature  (0.1  °F) 

4.  Air  orifice  size 

5.  Steam  orifice  size 

6.  Top  steam  header  pressure  (0.01  in.  Hg) 

7.  Bottom  steam  header  pressure  (0.01  in.  Hg) 

8.  Pressure  differential  across  the  steam  orifice  (0.01  in.  Hg 
under  H2O) 

9.  Core  upstream  pressure  (0.01  in.  H20) 

10.  Core  pressure  differential  (0.001  in.  H20  for  APc<  3",  and  0.01 
in.  H20  for  APc>  3") 

11.  Orifice  pressure  differential  (0.001  in.  H20  for  APo<  3",  and 
0.01  in.H20  for  APo  »3") 

12.  Orifice  upstream  pressure  (0.01  in.  H20) 

13.  Core  upstream  air  temperature  (0.01  millivolt) 

14.  Core  downstream  air  temperature  (0.01  millivolt) 

15.  Orifice  air  temperature  (0.01  millivolt) 

16.  Top  steam  header  temperature  (0.01  millivolt) 

17.  Bottom  steam  header  temperature  (0.01  millivolt) 

18.  Steam  orifice  temperature  (0.01  millivolt) 

19.  Condensate  weight  (0.1  lb) 

20.  Condensate  collection  time  (1.  sea) 

Record  items  6  through  18  three  times  in  the  same  sequence  and 
average  on  the  data  sheet. 

The  humidity  ratio,  H,  may  be  found  from  a  standard  psychometric 
chart  and  convert  temperatures  to   F   (0.1). 
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Figure  7.  Mass  rate  of  air  flow  vs. 

orifice  pressure  differential. 
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APPENDIX  II 
DIGITAL  COMPUTER  PROGRAM  FOR  DATA  REDUCTION 

The  program  is  designed  to  reduce  data  from  both  the  hot  core  test 
and  the  cold  core  test.   The  program  written  in  FORTRAN  60  and  is  shown 
on  pages  52  through  57  and  the  output  for  the  hot  core  test  is  shown  on 
page  58   and  that  of  the  cold  core  test  on  page  59  .   The  variables  used 
in  the  program  are  defined  in  the  program  glossary.   Variables  not  appear- 
ing in  the  glossary  are  internal  to  the  program  and  are  used  to  change 
dimension  or  to  define  special  groupings  to  simplify  programing. 

The  index  used  to  reduce  data  from  the  hot  core  test  or  the  cold 
core  test  is  given  by: 

Program  Index 
IHORC 

Hot  core  test  1 

Cold  core  test  2 

The  program  is  indexed  to  accept  the  various  fin  geometries  normally 

used  in  plate  fin  construction  and  the  indices   are  as  follows: 

Fin  Geometry  Index 
ISQTRIN 

Square  fin  1 

Triangular  fin  2 

Louvered  or  off-set  fin  3 

(NRe  =  infinity) 
•The  indices   used  for  the  fluid  metering  orifices  are: 

Air  Orifice  Air  orifice  Index 
BETA              d  (in.)  IB 
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.75  10.406  4 

.45  6.244  3 

.25  3.468  2 

.15  2.081  1 

Steam  Orifice  Steam  Orifice  Index 

BS  ds  (in.)  IBS 

.71  0.890  2 

.56  0.700  1 

The  indices  are  all  entered  as  fixed  point  constants.   The  number 
of  runs,  NOR,  and  the  run  number,  NR,  are  entered  as  fixed  point  con- 
stants.  All  fixed  point  constants  must  be  right  justified  in  their 
data  fields. 

The  constants  for  a  given  core  such  as  areas,  lengths,  etc.  are 
entered  as  input  data  only  once  and  are  the  first  three  data  cards 
beginning  with  AC  and  ending  with  DD.   The  fixed  point  constants  for  a 
particular  series  of  runs,  NOR,  IHORC,  ISQTRIN  are  entered  on  the  fourth 
data  card.   The  fifth  through  the  seventh  data  cards  are  for  run  one, 
eighth  through  tenth  for  run  two  and  so  on.   The  first  data  card  for  a 
run  contains  the  run  fixed  point  constants,  NR,  IB,  IBS.   The  second  and 
third  data  cards  for  a  run  contain  the  run  data  beginning  with  H  and 
ending  with  TC.   When  entering  data  for  the  cold  core  test,  only  the 
applicable  data  must  be  entered.   All  other  data  fields  may  be  left  blank. 

The  dimensions  of  input  data  must  be  those  given  in  the  program 
glossary  and  are  the  same  as  the  raw  data. 


46 


PROGRAM  GLOSSARY 

Program  Program  Nomenclature 

Symbol  Unit  Symbol 

AA  ft2  Aa 


c 


AC  ft2 

ALFA 

AFA  ft2  Afa 

AFR  ft2  Afr 

AFS  ft2  Afs 


fm 


AKFM  BTU/(hr  ft2  °F/ft)  K 

AS  ft2  Ag 

AWA  ft2  Awa 

AWS  ft2  A^ 

BETA  8 

BS  8S 

C  C 

cc  cc 

CO  c0 

COLJ  j 

C0S  Cos 

CP  BTU/(lbm  °F)  c 

P 

CRKD  K, 

ac 

CS  Cs 

DC  AC 

DCS  ACS 

DD  in .  D^ 
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DL 

DPC 

DPO 

DPS 

DSO 

ERROR 

FA 

FAS 

FB 

FBS 

FC 

FD 

FLA 

FLL  ft 

FLS  in. 

FKA 

FKS 

FM 

G 

H 

HA 

HG 

HS 

HS1 

HS18 

HS2 


in. 

*d 

in. 

H20 

APC 

in. 

H20 

AP 
o 

in. 

Hg  under  water 

APs 

in. 

ds 
Err 

FA 

FAs 
F 

Fs 
f 

fd 

in. 

*a 

^s 

BTU/(hr 

ft2QF/ft) 

kfa 

BTU/ (hr 

ft2oF/ft) 

kfs 
fm 

lbs/(hr 

ft*) 

G 

lb/lb 

H 

BTU/(hr 

ft2oF) 

ha 

BTU/ lb, 
m 

h 
c 

BTU/(hr 

ft2oF) 

hs 

BTU/ lb 

m 

hsl 

BTU/lbm 

hsl 

BTU/lbm 

hs2 
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IB 
IBS 
IHORC 
ISQTRIN 

NOR 

NR 

OD 

PI 

P2 

PB 

PO 

PR23 

PS1 

PS2 

PSO 

Q  AIR 

QSTM 

REB 

REFM 

REO 

RESO 

RH 

ROl 

R02 

ROAVE 


Index  for  air  orifice  diameter 

Index  for  steam  orifice  diameter 

Index  for  Hot  or  Cold  core  test 

Index  for  fin  geometry,  square,  triangular 
or  infinite  Reynolds  number 

Number  of  runs 

Run  number 

in.  d 


in.  H20 

in.  H20 

in.  Hg 

in.  H20 

in.  Hg 
in.  Hg 
psia 
BTU/hr 
BTU/hr 


ft 

lb;  m/f t3 
lb.  ^ft3 
lb  m/ft3 


N   2/3 
iNpr 

psl 

Ps2 

Pso 
Q  air 

Qsteam 

Nt 


*Reb 


N 


Ref 


N 


Red 


NResd 
rh 

Pi 
P2 
Pave 
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ROM  lb:  m/f t3 

ROO  lb;  p/ft3 

SDOT  lb   /hr 

pi' 

ST 

SQKD 

Tl  °F 

T2  °F 

TB  °F 

TC  min 

TFA  in. 

TFM  °F 

TFS  in 

TO  °F 

TS  °F 

TS1  °F 

TS2  °F 

TSO  °F 

TRKD 

TW  in. 

U  BTU/(hr  f?°F) 

UB  lb.^/Chr  ft) 

UC  lb  m/(hr  ft) 

UFM  lb  m/(hr  ft) 

UO  lbm/(hr  ft) 

USO  lb  m/(sec  ft) 


Pm 

Po 

s 

NSt 

Kd 
s 

•=1 

fc2 

cb 

time 

of  condensate 

collection 

*a 

tf 

5S 

^ 

ts 

tsl 

'■2 

tso 

Kdt 

I 

w 

u 

yUb 

yU   in 

core  for 

cold  core 

test 

Mt 
>to 
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vs  ft3/ibm  ys 

WC  lb  m  Weight  of  condensate 

collected 

WCH  lb  m/hr  wc 

WDOT  lb  bi/hr  m 

WK  BTU/(hr  ft2oF/ft)  kw 

X  x 

XC  X„ 

c 

SKD  Kd 

XM  Xm 

XMA  ft"1  ma 

XMS  ft"l  ms 

XNTU  NTU 

XS  xs 

Y  X 

YKE  Ke 

YNFA  nfa 

YNFS  Tifs 

YNO  n0 

YNS  ns 

YS  rs 

i&£  Kc 
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PROGRAM    SSHEAT 
OREAD    1,AC, AFR,AA,AWA, AFA, AS,AWS,AFS,FLL,TFS,TFA,TW, 
1FLA,FLS,FKA,FKS,WK,FD,DL,DD 

1  FORMATC8F10.0) 

READ  2,N0R, IHORC, ISQTRIN 

2  FORMAT(3I10) 

GO  TO  (3/5)* IHORC 

3  PRINT  4 


40FOKMAT(1H1 
1106H  RUN 


25X.32H 

MOOT 


HEAT 
G 


TRANSFER 
HUMID 


AND  FRICTION 
Tl     PI 


DATA,/// 


T2 


ERROR,// 
PSIA     F 
PERCENT,//) 


25H  ISOTHERMAL  FRICTION 
G       HUMID    Tl 
LB/HR  FT2  LB/LB     F 


DATA,/// 
PI     DPC 


2H  BTU/   COLJ    F  FACT    NRF      NR  B 

3107H      UB/HR   LB/HR  FT2  LB/LB     F 

4R  FT2  F 
GO  TO  7 
5  PRINT  6 
60FORMAT(1H1  25X, 

160H  RUN    MDOT 

246H       LB/HR 

7  READ  8,NR, IB, IBS 

8  FORMAT  (3110) 

OREAD  9,H,0D,DS0,BETA,BS,T0,T1,T2,TS1,TS2,TS0, 

1PB,PS1,PS2,DPS,P1,DPC,P0,DP0,WC,TC 

9  FORMAT  (8F10.0) 

PS1A=.490*(PB+PS1) 
PS2a=.490*(PB+PS2) 

PSAVE=(PSlA+PS2A)/2. 

TS=220.*(PSAVE-17.1)*2.66 

P2=P1*DPC 

T0R  =  T0«-459.7 

T0K=!0R*5./9. 

UO=,00  3527*TOK**1.5/(TOK+110,4) 

GO  TO  (20*21,22,23  ),IB 

20  CO=. 59446 
DC=.  009^5 
Cl=.5975 
OD=10.406 
BETA=.75 

GO  TO  25 

21  CO=. 59483 
DC=. 01037 
Cl=.5966 
0D=6.244 
BETA=. 45 
GO  TO  25 

22  CO=. 59868 
DC=.  01543 
C1=.6014 
0D=3.468 
BE1A=.25 
GO  TO  25 

23  CO=. 60480 
DC=. 05448 
Cl=.6128 
OD=2.081 
BE1A=,15 

25  C=C1 

B4=BETA**4 

POA=PB*,49Q-PO*,0361 

X=(DPO*,0361)/POA 

Y=l.-((( ,4l+.35*B4)*X)/1.4) 

XM=(1.*H)/(1.*1,607*H) 
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DPC 


PSI 


TS 


F  FACT    NR,// 


PSIA   PSI,//) 


R00=P0A*144.*XM/(53.3*T0R) 

IFCTO-68.  )12,12,11 

11  FA=l.+(T0-68. )*1.85E-5 
GO  TO  13 

12  FA=1. 

13  FB=1./(1.-B4>**,5 

DO  26  1=1,3 
WD0T=359.*C*FB*FA*Y*0D*0D*SQRTF(DP0*R00) 

RE0=l5,28   *WDOT/(UO*OD) 

C=CO*DC*SQRTF (10000. /REO) 
26  CONTINUE 

GO  TO  (29,40) , IHORC 
29  TB=(Tl+T2)/2. 

TFM=(TB+TS)/2. 

T8K=(T6+459.7)*5./9. 

TFMK=(TFM+459.7)*5./9. 

UB=.003527*TBK**1.5/(TBK+110.4) 

UFM=.003527*TFMK**1.5/(TFMK+110.4) 
AKFM=.001529*TFMK**.5/(1.*245.4/CTFMK*10.**(12./TFMK))) 

RH=AC*FLL/AA 

G=WDOT/AC 

REFM=4.*RH*G/UFM 

REB=4.*RH*G7UB 

XC=(1.+1.915*H)/(1.+H) 

pp=  XC*  241 

XNTU  =LOGF( (TS-T1)/(TS-T2)> 

U=CP*WDOT*XNTU/AA 

HS=2000, 

XMS=SQRTF(24.*HS/(FKS*TFS)) 

YNFS=(TANH(XMS*FLS/12.))/(XMS*FLS/12.) 

YNS=(AWS+YNFS*AFS)/AS 

HA  =  U 

DO  30  1=1*3 

XMA=SQRTF(24.*HA/(FKA*TFA) ) 
YNFA=(TANH(XM A* FLA/12. ) > / ( XMA*FL A/12 . ) 

YNO=(AWA+YNFA*AFA)/AA  i0  xXN. 

30  HA=l./(YN0*((l./0)-AA*(l./(YNS*AS*HS)+TW/(AWA*wK*12.)))) 

PR23=(UFM*CP/AKFM)**.667 

ST=HA/(G*CP) 

COLJ=  ST*PN23 

40  ALFA=  AC/AFR 

P1A=  PB*13.571-  PI 
P2A=  PB*13.571-  P2 
T1A=  Tl+459.7 
T2A=  T2+459.7 
R01=.097318*XM*P1A/T1A 
RO2=.097318*XM*P2A/T2A 

GO  TO  ( 41, 42) , IHORC  _  x  . 

41  R0M=.5*((TlA*R01+T2A*R02)/(TS+459.7-((T2A-TlA)/XNTU))) 

GO  TO  43  _ 

42  ROAVE  =  .097318*XM*(PlA+P2A)/(2,*TlA) 

ROM=ROAVE 

G=WDOT/AC 

RH=AC*FLL/AA 

TlK=(U  +  459.7)*5./9. 

UC=.0  0  3527*T1K**1.5/(T1K+110,4) 

REC=4.*RH*G/UC 
REFM=REC 

43  GO  TO  (50*45,44), ISQTRIN 

44  XKD=1. 
GO  TO  60 

45  IF(REFM-2000. >46,47,47 

53 


46 


47 


48 

50 
51 
52 
53 
54 

55 

56 
60 


TRKD=1.43 

GO  TO  48 

FM  =.049*REFM**(-.2) 

CRKD=1,090  68*4,*FM+.Q5884*SQRTF(4.*FM)*1. 

TRKD=1.+1.29*(CRKD-1. ) 

XKO=TRKD 

GO  TO  60 

IF(R£FM  -2000. )51,53,53 

SQKD  =  1.39 

GO  TO  56 

FM  =.049*RfcFM**(-.2> 

CRKD=1,09  0  68*4.*FM*.05884*SQRTF(4.*FM)+1. 

SQKD  =  l.+l.l7*(CRKD-l. ) 

XKD=SQKD 

CC  =.611+.045*ALFA  + . 344* ( ALF A**5 . 7 ) 

Y 

Z.. 

A=(ZKC-(1.+ALFA*ALFA))/R01 


70 
75 


76 


77 


80 


KE=1.-2.*XKD*ALFA+ALFA**2. 
KC=(l.-2.*CO(CC*CC)*((2.*XKD)-l.))/(CC*CC) 

A=(ZKC-(1.+ALFA*ALFA))/R01 

B=(YKE+1.+ALFA*ALFA*(1.+(4.*FD*DL/DD) ) )/R02 
FC=   (RH*R0M/FLL)*((4.3255E+9*DPC/(G*G) )-A-B 

r...  a  a   -   niAj.  r,  1  £  * 


104)  , IHORC 
76), IBS 


P1AA  =  P1A*.0361 

DPCA=DPC*,0361 

GO  TO  (70, 

GO  TO  (75, 

COS=.6090 

DCS= .0410 

DS0=.89 

8S=.71 

GO  TO  77 

COS=.6085 

DCS=.0250 

DS0=.7Q 

BS=.56 

CS=.61 

BS4=BS**4. 

FBS=1./(1.-BS4)**.5 

FAS=1.003 

DPS0=.454*DPS  0 

PSO  =  17.1+(TSO-220.  **.376 

XS=DPSO/PSO 

YS=1.-(.41+.35*BS4)*(XS/1.3) 

DPSW=DPSO/,0361 

VS=24.01-(TSO-218. )*.4 

USO=8.25E-6+(TSO-200.)*1.45E-8 

DO  80  IS=1,3 
SD0T=359.*CS*FSB*FAS*YS*DS0*DS0*SQRTF(DPSW/VS) 

RESO=.004244*SDOT/(USO*DSO) 
CS=COS+DCS*SQRTF(10  0  0  0./RESO) 

CONTINUE 

HSl8=ll58,0+(TSl-230.)*.5 
HS1=HS18-(PS1-18.0)*.35 
HS2=1155.3+(PS2A-19. )*( .934) 
HC=193.4*(PS2A-19.)*2.64 

WCH=WC*60./TC 
QAIR=WD0T*CP*(T2-T1) 
QSTM=WCH*(HS1-HC)*SD0T*(HS1-HS2) 
ERROR=(QAIR-QSTM)*100./QAIR 
99  PRINT  10  0,NR,WDOT,G»H,T1,P1AA,T2,DPCA,TS,HA,COU#FC,REFM,REB, 
10Q0FORMAT(I3,2X,F6.1,2X,F7.1,2X,F5.4,2X,F4.1,2X,F5.2,2X,F5.1,2X, 
12x,F5.1,2X,F4.1,3x,F6.5,2x,F6.5,2x,F6.1,2X,F6.1,2X,F5.1,//) 

GO  TO  110 
104  PRINT  105,NR,WDOT,G,H,T1,P1AA,DPCA,FC,REFM 


ERROR 
F4.2, 


1050FORMATII3.2X.F6.1,2X,F7.1.2X.F5.4,2X,F4.1.2X,F5.2,2X,F4,2.2X,F6.5. 

12X,F6.1,//> 
110    IF(NR-N0R)7*112,112 
112    CONTINUE 
END 
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APPENDIX  III 
GRAPHS  FOR  APPROXIMATE  EQUATIONS 


Graphs  for  the  approximate  equations 
used  in  the  experimental  method  and 
equations  section  are  given  on  the 
following  pages. 
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Figure  8.   Thermal  expansion  faotor 
vs.  temperature. 


250 


61 


1.00 


0.95 


0.90 


0.85 


0.80 


0.75 


0.70 


0.65 


0.6o 


Figure  9.  Contraction  coefficient  vs.tf" 
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Figure  10,  Saturated  steam  pressure 
vs.  temperature* 
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Figure  11.  Saturated  steam  dynamic 

viscosity  vs.  temperature. 
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Figure  12,  Saturated  steam  specific 
volume  vs.  temperature* 
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Figure  15,  Saturated  steam  enthalpy 
vs.  pressure* 
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Figuro  14.  Superheated  steam  enthalpy 
vs.  temperature. 
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Figure  15.  Saturated  liquid  enthalpy 
vs.  pressure* 
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